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ABSTRACT 


An experimental investigation has been carried out 
to demonstrate that the well known sheath/convection and 
sheath/diffusion theories are not applicable to measure- 
MNentSeOt eaelegclectrongdens 2 ysO act tues Gh res sine 
flowing plasmasusing an ion probe. 

It has been found that at low ionization densities, 
the current is no longer space charge limited but is in 
facta limited by theselectric tiela generated by the prope. 
Experimental measurements with a cylindrical probe at 
densities down to Mee ae show good agreement with the 


Concept.oLea 1.26.10 limited current, 
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Tole OFeoYMBOLS 


Ry = Electric Reynolds number 
NaS Plasma convection velocity 
= = Probe radius 
D; = Positive ion diffusivity 
eQXt 
n ae 
ae Debye Length _ € 

Probe Radius rs 
oe Permittivity of free space 
k = Boltzmann's constant 
oe = Electron temperature. 
a = Electron density 
e = Electronic charge 
xX = ev_ 

kT, 
V = Probe potential (Negative) 
eee LOmice MOM LL Ly: 
E = Electric field 
the ee Ion density 
Q = Charge per unit length on the probe 
Vv == Velocity Of tons relative to) the probe 
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a = Probe current per unit length 


d = The angle between the velocity and positive vectors. 
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CHAPTER. I 


INTRODUCTION 


The development of accurate, well understood plasma 
diagnostic techniques (i.e. determination of plasma char- 
acteristics, mainly electron and ion density, electron 
temperature, recombination rate etc.) is essential in plas- 
ma research. Any progress in the understanding of natur- 
ally occurring plasmas or in utilizing existing techniques 
in the production of plasma will put this science another 
step ahead. Primary sources of naturally occurring plasma 


include the ionosphere and flame plasmas. 


1.1 FLAME AS A PLASMA SOURCE 


The convenience of the production of flames makes them 
a very important source of plasma. They have a temperature 
range from 1000°K to 4000°K and a pressure range of typi- 
Callvyet rom jl) torn to lO0 atmospheres es Ine electron concen— 


tration of a flame can vary from 10> fm? to aye yee with 


6 4 


to 3 X 10 ‘m. 


oud Ore 


Debye shielding lengths of the order of 3 X 10= 
The usual characteristic flame dimensions of ie 
are sufficient to allow the flame to be treated as a true 
plasma as. distinct from a simple ensemble of ions and elec- 
trons. The collision frequency of electrons and neutral 


i -13 
molecules in flame gases is of the order of v = 10 1 IM]/sec 


3 
where [M] = Concentration of neutral molecules per meter’. 
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In an atmospheric pressure flame v ~3 x Ogg /secnle 
Because of these qualities research workers have used 


flames as a convenient source of plasma for many experi- 


ments. 


Various techniques have so far been proposed for the 


study of flame plasmas namely, 


thes Microwaves 
Dr Electrostatic Probes 
ae Mass Spectrometry 


4. Spectrophotometry 

ae Electron Cyclotron Resonance 

6. Radio vrrequency Coils 
The first three methods have met with a good deal of suc- 
cess. 


q) suffers from a drawback be- 


The microwaye metho 
cause it needs space to place the flame within the cavity. 
The spatial resolution is thus reduced by the presence 
of large holes in the cavity walls. The holes in the cay- 
ity which are necessary for the entrance and exit of the 


hot flame gases also permit the electric and magnetic 


fields to leak out of the cavity. 


The limitations of the mass spectrometer are associated 
with the problem of accurate sampling of the flame gases. 


Although much work has been done to alleviate the problem, 
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it is at its worst at high gas pressures and extreme care 


must be taken in identifying true flame ions. 


1.2 ELECTROSTATIC PROBE AS A PLASMA DIAGNOSTIC 


Electrostatic probes are one of the fundamental tech- 
niques in plasma diagnostics. In 1924 Langmuir developed 
this technique. The probe is a metallic wire inserted in- 
to the plasma and biased positive or negative with respect 
to the plasma with the help of a power supply or battery. 
The current-voltage characteristics provide information 
about the condition of the plasma. For a wide range of 
conditions the presence of the probe does not perturb the 
plasma, Thus the probe does not have any effect on the 
quantities being measured. The electrostatic probe offers 
good spacial resolution oyer a wide range of conditions, 
Although the probe is very easy to construct, the inter- 
pretation of the results is not so simple.” The dzritireulty 
arises because probes are boundaries to the plasma. The 
governing equations of motion of the charged particles in 
the plasma do not hold in the boundary region. The con- 
dition of quasineutrality which 2S true “for the plasma is 
not valid near the boundary (which is known as a "sheath") 


where the electron and ion density differ. 


Langmuir and H.M. Mott smith ©) assumed a thin bound- 


ary layer around the probe and noted that the quasineutral 
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equation could be used up to the 'sheath edge! which in 


practice had a well-defined position. 


In recent years scientists have tried to understand 
the behaviour of a probe in a high pressure plasma. Schultz 
(4) 


and Brown tried to interpret the positive ion saturation 


region of a probe at low pressure, taking into account the 
directed current at the sheath edge, Good agreement was 
found between the theoretical and experimental results 

at low pressure. At high pressure the probe collects a- 
bout twice the current predicted by the theory. Their 
theory was based on the assumption of a sheath region 
around the probe. Su and tam ©? have presented qa contin- 
uum. theory fora negatively biased probe in a slightly 
ionized high neutral density gas. They assumed a slightly 
ionized plasma, which was quiescent. They have given two 
explicit forms of current-voltage characteristics, one for 
yery negative probes and the other for probes at close to 
plasma potential, Both of these are based on the assump- 
tion that the probe radius is large compared to the Debye 
length. Soundy and Williams ‘2) developed a probing tech- 
nique along the lines suggested by the theory Ole oU.ald 
Lam, but the current collected bythe probe was 100 times 
greater than the theoretically predicted current, A nunp-= 


ber of theories have been put forward to explain this dis- 
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1.3 ELECTROSTATIC PROBE FOR THE DIAGNOSTICS OF FLOWING 
PLASMAS 


It has been found experimentally that the electro- 
static probe theory works satisfactorily in static plasmas. 
However, the theory does not work in most atmospheric 
plasmas since in fact such plasmas must usually be con- 
sidered as moving plasmas, because of their inherent mo- 
tion. There has been considerable interest in the appli- 


cation of probe theory in this type of moving plasma ‘©712) | 


(6) 


Lam first Looked anto the caseyot nighypressure 
moving plasmas For a probe, the condition of ion collec- 
tion reduced-to-that considered by Su and Lam for the case 


of electric Reynolds number Rz<< 1, where the electron and 


ion temperature are approximately the same. 


* =1 
EI Bias Be, 
where 
es plasma convection velocity 
Xp = probe radius 
Be = positive ion diffusivity 


This theory is only applicable under the following con- 


ditions: 


1. The mean free path of the charged particles is 
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much smaller than the thickness of the sheath. 
2. The Debye length is much smaller than the thick- 


ness of the boundary layer adjacent to the body 


surface, 


The prediction of Lam was that the probe acts as a 
Sink to the charged particles which lose their charge by 
recombining on the surface of the probe. The electrons 
will hit the probe more often than the heavier ions due 
to their large thermal velocity. Immediately next to the 
body surface the number density of electrons becomes too 
low to carry the ions and the ions' own diffusion motion 
takes over. Outside this ambipolar diffusion region, con- 
vVect1OnM eLtects are dominant, and these ectronpanidelonmacens 


Sities are uniform and equal. 


C3) 


De Boer and Johnson have derived the sheath thick- 
ness and current collected by the probe which depends on 
the conservation of ions entering the sheath from the free 
stream. They haye assumed that the boundary layer between 
the quasineutral region and the sheath edge is very small. 
Aestrong electric field acts at the ion sheath and thevefq 
fect of the electric field on the ion flow is very small. 
The free stream electric field cannot be precisely zero, 


since some field strength is required to drive the elec- 


tron towards the positive electrode. The ions having been 
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convected through the sheath edge are assumed to change 
direction abruptly and then move directly to the probe, 
under the influence of strong electric fields, in the 
sheath. Again since the plasma is flowing, a hydrodyna- 
mic boundary layer is formed around the probe. If the 
electric sheath is much thicker than this hydrodynamic 
boundary layer, convection within the sheath may be ne- 


glectea‘1*) . 


The mechanism of ion flow within the sheath 
can be treated as one dimensional, which greatly simpli- 
fies the task of solving the sheath equation. At low 

ion density the sheath becomes thick compared to the hy- 
drodynamic boundary layer. For sufficiently low ion den- 
sity, the sheath can fill a hollow probe, which then 
Servesmacedutotal collector sor al litonseconvectedminros.tes 
entrance area. In describing the motion of ions ina 


ee concluded that convection simply adds 


sheath, De Boer 
a constant horizontal component to the ion motion in the 
major part of the sheath, The predominant electric ef- 
fect is in the vertical direction and not in the horizontal 


direction, but the electron density from probe measure- 


ment was 50% more in De Boer measurement. 


Clements and Smy divided the high pressure moying 
Zee 
probe theory into two regions namely, Ro So ang 
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e 


os 


7 ai ie 7 . 
weit pes eae 


eT 7 a en 
wih panes? 


Saari * wet alot 


Se i ee) oer, 
Dsdedd ae re lt! te ba fea aa hee 
sifened ydainsp. alte Spabiis Anity ORG ote he 
wel-gA cea dada dies iit Ses-oe kets ey ie 
“an 8° 2 Ajeru> & iti Sui sx saat i aes ae 
ae Ch Wel eile. HA ai fratinvia ro ae 
fe RL wheg apie pie “er eucehs ane “i. 
O54 Geni bby sce; owh- Akh ty guess Sa thSoas “a 4 ‘wwe a 
Pw A ee 
ee ee ee ee * 
9 wale FEA AD, 93, 9 gr, deta rshe ee 
sie ste Pails Sasabed ss Hp ho ian kel io oteq SOR | 
(eemeetac’ ef 6! TH PAs healer tLe Lesey? et a4 <n | a 
er eruluit’ ere: om Aleits ee ad? ua pol 
Sretemdsaan 41267 Ba kt ow Oye ‘eon 1m? 


petvon weeny Mat GA? ey 4s et ep ssi hD 

alee arte pad oy ru vit a 

ibors iambbe~ a jail 6 - 
' 7 


7 7 


- - 
a 
a 


Cako 
2 vligiaen ©) 

ner Z 

€ p 
Coe Permittivity of free space 
k = Boltzman's constant 
ae = Electron temperature 
n, = Electron density 
e — Eh Lectronic charge 
i = Probe radius 
x = ev/kT . 
V = Probe potential (negative) 
ies Lone mobi Uunty: 


For Ryo x7 >>1 Clements and Smy haye shown that the 


supply of ions to a negatively biased probe arises pre- 
dominantly from hydrodynamic convection of ions into the 
probe sheath. If the sheath remains thin compared to the 
Peobe ecadlus;.the’theovreticalscurrenteperaunie Jengthagco 
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I = 5.3 (n,ev_) . Zale 
For R, a? <<1, the supply of 10ns to the probes segovern— 
ed by diffusion of ions through the electric boundary lay- 
er surrounding the probe. The theoretical diffusive cur- 


rent per unit length to a cylindrical probe was derived 
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By following the guide lines presented by Lam, Cle- 
ments and Smy ee came up with a simple expression for 
the transition between diffusive and sheath convective be- 
haviour in the absence of recombination. They found this 
transition to be a reasonably broad one, which in fact 
covers a range of about two orders of magnitude of the 
variable (Re Ao yh Additional complications arise 
from the cold dense gases in the hydrodynamic boundary 
layer. When this boundary layer is thick in comparision 
with the sheath, then the boundary layer resistance has 


a large influence on the probe current. 


1.4 THE NATURE OF THE PRESENT WORK 


There#is no definite theory to calculate the electron 
density from the saturation ion current under a variety 
of plasma conditions. The dominant supply of ions to a 
negative probe arises from their convection or diffusion 
into the sheath surrounding the probe. In the simple ™ 
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model proposed by De Boer and Johnson Kulgein 


Clements and smy 649) the ions had been convected through 
the sheath edge. Once they cross the sheath edge, the 
ions change direction quickly under the Lnoluence -OLe the 
high aeseredie field which occurs in the sheath and move 


towards the probe. This theory has been applied by others 
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PosdEEtecrent= probesa (cylindrical, Spherical, and stagna- 
tion probe) and sheath configurations (planar, cylindri- 
cal). Good agreement has been found with existing exper- 
imental data for ionization densities in the range of Oe 
LOE The sheath convection theory gives good agree- 
ment, but a significant error arises in the case of sheath- 
diffusion theory (De Boer (13) vl10x, Clements and smy (+8) 
vl00x). This particular mode of operation, the so called 
sheath/convection mode, is particulary desirable because 
for many geometries the theory is simple; its assumptions 
are reasonably well satisfied and the measurement is per- 
formed on an ‘undisturbed plasma!, since the sheath edge 
which abstracts ions from the plasma. extends beyond the 
hydrodynamic and thermal boundary layers surrounding the 
probe. As has been pointed out above, the application of 
the sheath/convection model is limited by yarious para- 
meters. This thesis describes an experimental investiga- 
tion of probe behaviour at ionization densities down to 
ae ie in a region where the sheath/convection model 
loses validity and therefore another model must be used, 
The theory developed for this new region of ultra-low 
ionization densities (for laboratory plasmas) also is 


shown to have application to the measurement of ion den- 


sities within thick sheaths or clouds of positive ions. 
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Hout (19) and Sonin 29) have developed a theory for 


the blunt nosed probe which has some similarity with the 
present theory. Their theory does not consider cylindri- 
cal or spherical geometry and takes no account of the 

sheath effects. The experimental measurement shows good 


agreement with the theory presented here. 
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CHAPTER LL 


THEORY 
2.1 THE LIMIT OF SPACE CHARGE AND FIELD LIMITED CURRENT 


The simple sheath/convection model by De Boer and 


siciee (21) 


, and Clements and Smy TS Shown inetd. ele 
In the case of a cylindrical probe in atmospheric plasmas, 
the sheath radius is much greater than the probe radius. 
The essential point of the sheath/convection model is that 
all ions convected into the sheath eventually arrive at 
the probe due to the strong electric field in the sheath. 
The field E = 0 at the sheath edge (the boundary of the 
plasma electrons), and as soon as the ions enter into the 
sheath their motions are determined by the sheath elect- 
ric field. The sheath radius then adjusts to such a value 
(x0) (assuming a cylindrical symmetry) that the’ current 


to the probe is space charge controlled. The ion cur=- 


rent/unit: length at a distance © is given by 
rhe od Cee le nu; Be 15); 


where ye ion mobility, & =) there lectricet tel soteche 


electronic charge, De ionization density. 


probe, e 


From Poissons equation, in cylindrical coordinates 
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Schultz and Brown have integrated this equation ap- 


plying the boundary condition E=V=0 at r = ie and V = 


nS acer = XD and i= I giving rise to the expression 
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where a Permittivity of free space, 


It is noted that in deriving the asymptotic expres- 
Sion above, a factor of 1.37 in the 'log' term has been 
neqglectediforisimplicityslaewsupply Of ftonsmtosEehegoucer 
edge of the sheath will now be generated by convection 


LOD Tot and Be >l and will consequently be given by 


= 4) 
ie 2ron.eVe (4) 


R_ = electric Reynold's number = AN EE yee 


= flow velocity of the plasma, 


We 
k = Boltzman's constant 
T. = electron temperature. 


It is apparent from equations (3) and (4) that the sheath 
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thickness lo will adjust so that the linear conyection of 
ions into the sheath is just balanced by the radial space 


Charge limited flow of ions into the probe. Eliminating 


ry from equations (3) and (4) Clements and smy (24) ob- 
tained 
aie Sug ee tee 
V = log ( ————————_) 
a . 
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At low ionization density or at high flow velocity 
it is observed that the fields within the sheath may not 
be strong enough to ensure that a large portion of the 
ions entering the sheath do indeed reach the probe; many, 
in fact, will be swept on by the flow back into the plas- 
Mase hiLcmsituatLon=seindicated invtig.  (l-b)iee nem ons 
between AAswill only contribute to the probe current, sthe 
others will not contribute to the probe current although 
they remain in the sheath for a period of time, It can 


be decided which of these regions best describes a giyen 
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probe situation by considering the situations (1—b),.1in 
this case the motion of most ions within the sheath is 
largely unaffected by the sheath electric field. Con- 
sequently, the ion density remains close to the value 
in the plasma and with a cylindrical sheath of radius 


"ry! with potential 'V!' we have 


DiGi tan ree ne where r_>>r 
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(6) 
The sheath electric field then can be calculated and is 


given by 


O Ae (7) 
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FIGURE (l-a) Sheath/convection model from ref-22. All ions convected 


into the sheath eventually arrive at the probe. 
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FIGURE (1-b) Field limited model. Only the ions in the band AA reach 
the probe, 


Meee sn) planar sheath model. Eo is the local field. 
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It is noted that the assumption of the ion motion 
will hold if the product of this electric field and the 


ion mobility « ion convection velocity. 


y,E << Ve 
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Hs = <<V 
1 4 is £ (8) 


Converting the equation (8) into non-dimensional para- 
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2.2 THIN PROBE CYLINDRICAL GEOMETRY 

If the motion of most of the ions within the sheath 
is not affected very much by the sheath electric field, 
then the ionization density everywhere in the sheath 
apart from the region close to the probe will be uni- 


form and equal to the ion density in the plasma. 


Let Q be the charge per unit length on the probe, 
which is equal and opposite to the total charge of the 
ions in the sheath. The electric field in the region 


close to the probe will be given by the expression 


(10) 
where r is the unit radius vector to the probe. This 
region is the most important region, since it is only 


from this region that ions can be removed from the gas. 


If an ion moves with velocity vector v relative 
to the probe and the velocity of gas relative to the 
probe is Vi¢ then the velocity of the ion relative to 
Elemis al seGya— Ve) , and which has been produced 


by an electric field 
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From equation (10) we have 


Thus the path of an ion relative to the probe is des- 


cribed by the velocity vector 
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Wencan convert the equation (%2) in termsgory the non- 
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dimensional quantities using X = <n L and 
y2me Ve - 
Y = ———— ’ where x and y are the co-ordinates parallel 
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and perpendicular to v, respectively, with the origin at 


te 
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Thus the component of v in the y-direction, from figs (l=): 
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The component of V Simethe x-direction, from fig (1—d); 
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Equation (15) has a solution given by, 


is 


Y - tan = Constant (16) 


IK 
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Plots of this equation appropriate to three dif- 
ferent values of the constant) are Shown in) tid+(2)-s Tne 
solid line shows the trajectories of ions whose initial 
value at x = - © is Yi= one The ions whose initial value 
at x = -~ is less than he will finally meet the probe. 
The ions whose initial values are greater than oe will 
pass on towards x = += , Inspection of equation (16) 


above shows that this situation occurs when the constant 


= 0 
Y = fan > = reboyel yeh “Ne MP It] y.€ eo th 
alte its sO ¥ shea. (77) 
O fe) 
Only those ions having the value of Ye Y> aha 
at x = -® will reach the probe, and thus the ion cur- 


rent /unic length of probe within the band Y= a will 
be 
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Asymptote to Solid Line=1] 


— eet 7 The Boundary Conditions 
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RPIGURE (2) “The trajectories for ions having different” impact 


parameters solved from equation (15). 
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(20) 


sibhe Che oe te pee ea 1 
dMiexefores the ion current/unit lengthiis 
n,eQu; 


a es ace | (21) 
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To obtain a current-voltage relation for such a 
probe it remains to relate the charge Q to the voltage 
applied to the probe. In the basic model we assumed 
that most of the ions within the sheath were essentially 
unaffected by the probe. Thus, for the purposes of 
calculating the probe voltage the sheath can be consid- 
ered as a cylinder of radius ro uniformly filled with a 
charge density of noe with no electric field or potential 


at its surface. 


The negative charge on the probe will be equal and 
opposite to the positive charge in the sheath. We see 
thatthe charge Q on the probe will be = -1r,’n.e which 
together with the volume charge of the sheath will 


produce an electric field at any POln tie 
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the sheath is 
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Pieyeamebrackets in equation. (22) can be approximated 
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But from equation (21) we get 
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Putting the value of Q from equation (23) one gets 


en _u.. LO ee VASO. 
I. = = = men, x cs ( = >) 
fe) BP Ses 
2 
ZL5enaG 
= ren iup.V 2 Mose 
Oat EGV 
2) s3ener 
= 2 nren_y.V = 0.3 
Oa enV. 
fe) 
Again, from equation (22) we get 
2 2 2 2 
Ne ry Lo Coie 5 
eg ec ee 
e ”) A ty 7, 
O r p is 
Pp 
ner 4 he ae 2 
= = O O {log pe) - P| 
2€ 1g e) 
Oo ae 
But QO == TX, ne 
2 
Q Tpek yy 
V= CEO mag ieee at 
ZTE, 2 Pig 


Spee mies tee 
e) Pp 


26 


(24) 


sn 

; : . 

ihe 

— 
ad =e a 


2 P | = 
qv ~ 
la 2 
— "4 ‘ ~~ . oa i 
> al bial a S ‘els ‘ 
> od a t 
_ rm On t : 
ty & 
& ane a = 
a ee 4 7 ee ; - 
7 ois 
* + “ ars - 
" ee 
ay = 
7 | = : —_ 


eae 


2Te0V 
(Oy r _ 
ey > 
arose 1G 2 
p 
Putting the value of Q in equation (21) one gets 
: 27en ou; 
p 18 
Oh atl (26) 
log ae > 
p 


2.3 THICK PROBE CYLINDRICAL GEOMETRY 


With high flow velocities or with large diameter 


probes the distance ve becomes very small compared) to 


Ene wprobpe radius. <«<In\this *case, flow is dutferent., | tne 
electric field, equal to oe , Will be constant while 
O 


the flow velocity considering the inviscid flow varies 


with position and is given by 


Ma Ces sing C2) 


where ¢ is the angle between the velocity vector Ve and 
the radius to the point concerned as shown in fig (l-e). 


The ion moves through an angle df? in time dt. 
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and during this time the ion will move radially inwards 
under the action of an electric field by an amount dr, 


therefore, 


de Quy 
dt ~ Qe lr 
B 
ut dr _ dr dt Q Ws 
do dt do Ate Ve sing 
Qu. 
Be ~ Gte_v ant fo 
OL 


(15) in a Similar situa- 


Following Clements and Smy 
tion involving the sheath convection model, it is neces- 
sary to modify equation (25) in the flow and an extra term 


is to be added for the inherent radial flow velocity. As 


a result equation (29) becomes 
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which has the solution 
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The current/unit length to the probe will be 
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- 4 = = Ouvenre 
. 2 (2v, r sind yao) eae Spe kOe (31) 


2.4 SPHERICAL PROBE GEOMETRY 


ine thes case™orea spherical probes the mnviseldenlow 


velocity which varies with position is given by, 


v= 3 Wes Woh 
oy hig ean) 


In time dt, the ion will move through an angle dd and 
during this time will move radially inwards under the 
action Of an e€lectric field. Considering the inherent 


radial Elow velocity in this Situation, one gets, 


dr Qu; 2. 2r cosd¢ 
do 4 Bye cine 6 sin¢ (32) 
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The current to the probe is given by, 
te = nev2msing i ee 
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But, one knows that 
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Bucting the value of 0 in equation. (34) one gets 


I = ae Rol ea’ (36) 


2.5 PROBE MEASUREMENT IN AN EXTERNALLY GENERATED THICK 
SHEATH 

The motivation for these calculations arose froma 
somewhat different experimental situation, It is possible 
to determine the ion density profile by means of measure- 
ments with a small probe in a thick planar sheath. This 
SLutuationels shown in Gig. G=—c).. Im this tcype On sicath 
region there are only neutral gas molecules and ions in 
an electric field. By considering the probe and its bias 
below its surroundings to be small, the sheath electric 
field can be considered uniform over the region of influ- 


ence of the probe and the ion density typified by a value 


Oe 


The total electric field on an ion is given by 
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which causes an ion drift velocity of 


2 Ors 
Vieaar it ontabintems md (37) 
The relation between Q and V for the static sheath 
is, however, different from that for the flowing plasma. 
It is assumed that the probe does not affect the boun- 
daries of the sheath and so the usual in 'vacuo! expres- 


sion applies, 
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where 2X5 veo the length, typical of the dimensions or ethe 


sheath. Putting the value of Q in the equation 
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which is similar to equation (26). The expression for a 
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spherical probe in a static sheath is the same as equa- 


eon m3 6)e 


In deriving equation (39) we have assumed that the 
ion current near the probe does not distort the electric 
field. At high densities such distortion will occur and 
the current to the probe will be limited by the space 
Charge. Effectively the voltage applied between the probe 
and its surroundings will drop off across the space charge 
region and only the ions swept into this region will reach 
the probe. In this situation equation (5) will apply be- 


cause of the sheath convection situation. 


In this section four expressions for the current to 
cylindrical and spherical probe under different conditions 
are derived. These formulas can be applied at lower den- 
sities in flowing plasma than before and also can be used 
to carry out direct measurements of ion density in sheaths. 
The most important differences between the conventional 
sheath/convection expressions and these derived expres-—: 


sions are the absence of any dependence on flow velocity. 


Physically, the probe can be considered as having 
two limits on its influence on the plasma. The first (ap- 
propriate to the electric field that would exist in vacuo) 


is dictated by the ratio of drag forces to electric forces, 
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and the second corresponds to the probe length which 
determines the actual probe currents. Spi ibllkeweiie she 

is the smallest of the field limited or sheath/convec- 
tion current which is appropriate. Comparing the two 
Ee pGesc lonsmtorm the cylindrical current7equation (5) 
(thick sheath) and equation (26) (field limited),we 


observe that the ratio 
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and so we might expect from the discussion leading to 
Ro 
would define an approximate 


relationships (9)) that 
x 
boundary between the field limited and sheath/convection 


regimes. 
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CHAPTER ITI 


EXPERIMENT 


TN TRODUCTION': 


An experimental investigation has been carried 
outs to prove that at low ionization densities the 
probe current is no longer space charge limited but 
is, in fact, limited by the electric field generated 
by the probe. A Meker type of burner was used to pro- 
duce a high pressure flowing plasma by burning pro- 
pane and air. Potassium, an alkali material, was used 
to increase the ionization level in the flame. The 
current, which gave information about the condition 
of the plasmawas collected through an air cooled 
brass probe. An air cooled gridwas used to create 
a planar sheath in the plasma. The flow velocity 
and the temperature of the flame:wre measured to 
Scaremthe mobulityeoteions. *The details ofsthe tex= 
periments and their setup are described in different 


sections of this chapter. 


31% THE BURNER SYSTEM 


The flame is the source of plasma in the exper- 
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iment. The Meker burner, which produces a uniform and 
Stable plasma, has a small nozzleshaped orifice of 1.5 
mm diameter for the gas with an air inlet located just 
above this orifice. The air is mixed with the gas in 

a long burner tube before combustion. The long burner 
tube reduces the turbulence due to the gas jet that is 
flowing through the burner. The optimal ratio of length 
to diameter of the burner “ree be 6:1, whereas the 
ratio used in the experimental burner is 5:1. The Meker 
burner used in our experiment gives improved air flow 
and so achieves a higher temperature and better ion- 
ization aes The metal grid at the top of the burner 
prevents flashback and thus achieves a more stable flame. 
The size of the holes in the grid is not very critical, 
although a large hole size may increase the risk of 
flashback. Lf the holes are too small and close toyeach 
other, there is a tendency for the normally separated 
flames to fuse together and lift off the grid surface 

to give a conventional bunsen shaped cone. The central 
apertures sin the orid give slightly different combustion 
conditions, due to the lack of secondary air reaching 


that region. The grid also acts as a heat sink. Grid 
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Supported flames show much better stability <2) owing to 
the damping effect which the screen sink has on the 


front movement of the flame, 


Although the front movement of the flames is re- 
duced=by.the grid Tt is ‘still the major problem in the 
burner aH. The flame also moves about due to the 
poor flow regulation. The problem of oscillation of 
the flame is reduced to a great extent by the grid. 
The vibrationaly induced eddies can be eliminated by 
enrixed@platform.» |The instability offa rising not) gas 
column is responsible for the flickering in the tip of 
the flame. It is best alleviated by introducing a 
chimney which captures the hot gas column before the 
instability develops. This provides a reference hot 


boundary and adds_to the stability of theltlame. 


The experimental setup of the burner is shown in 
fC) ae Lhe iign pressure propane passed through 
a regulating bulb to a propane tank of 3" diameter and 
12" length, The propane enteredthrough the bottom 
inlet and came out from the top. The tank helped 
keeping a continuous and steady flow Of Propane CoO. the 
burner. The gas then enters « the burner through a 
small orifice at the bottom and wasregulated with a flow 


meter placed in between the tank and the burner. Air is 
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also passed to the burner through two pressure regulating 
valves, an air tank and a flow meter. ‘Two chimneys are 
placed at the top and at the bottom parts of the flame 
using clamps. It was found experimentally that the 
chimney at the bottom is the most effective in keeping 
the flame stable. Photographs of the flame without 

a chimney, with a chimney at the top and also at the 
botton@arcyshown ins figs. (4) and) (5). ethe function of 
the exhause pump at the top is to pump out the burned 


gases. 
3.2 FLAME IONIZATION 


It has been found that the most important charged 
species present in flames are electrons and corresponding 
amounts of}jpositive ions. There are also small traces 
of negative ions in meas oO The concentration of 
ions found in the reaction zones of unseeded premixed 


hydrocarbon/air flames at pressure between 2 to 760 
15 918 mo (25728) | 


mo Ho Wiessin = the range yor 70 a Ak ve The 
+. 
only ions common to all hydrocarbon flames are CH, , ¢ ts 
Supe <a 
CHO’, and H0. Ofe these, H,0 is found later than the 
(22) 


others and therefore is unlikely to be a primary ion 
In the case of the premixed flame, the rate of ion 
generation increases steeply with final flame temperature 
(typically by 200% per 100° C). The rate of ion 
generation also varies with the air/fuel ratio and 


is a maximum at the stiochiometric ratio. The ion 
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concentration of a:propane/air flame is found to varyawith 
pressure and is a maximum at atmospheric pressure (2°) , 
The thickness of the reaction zone of the flame depends 


on pressure and is a minimum at atmospheric pressure (79) , 


Clements and smy ‘18) gouna that in a propane/air flame, 
the ionization is maximum at the reaction zone and 
falls with distance from the reaction zone. The ex- 
perimental measurement shown in fig.(6) agrees with the 


results of Clements and smy (18) | 


The experimental 
measurements were carried out in the range of 2 cms 
to 10 cms from the top of the burner and it was observed 


that the ionization remains fairly constant in that 


region, 
Se LONLZATION IN THE YP RESENCH OF ADDITIVES 


The ionization level of flames can be increased 
by adding small traces of substances having low ion- 
ization potential. The rate of ionization in flames 
is found to be a decreasing function of ionization 


potential Li<Na<K< Rb<Cs. 


The positiye ions produced (together with electrons} 
tend to consist of simple molecular compounds, and are 
often monoatomic metal ions. Consequently, the corres~ 
ponding thermal data (ionization potentials) are known. 

It is also easy to add nebulized solutions of metal 


salt to the gas supply and to vary these over wide 


pace 7 . oy 
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ranges by changing the concentration of the solution. The 
ionization of the metal is not restricted to the re- 
action zone, but occurs throughout the burnt gases. In 


the reaction zone, the dominant process of ionization 


is “chemi ionization". The dominant arose ion see 
is produced due to chemi ionization, which reacts with 
the metal atoms in a process shown below (24) 

A+ H,0%__»at 44,0 +H 


Where A = alkali metal 

This transferred chemi-ionization is capable of producing 
values of in] well above those appropriate to thermal 
equilibrium which persists into the reaction zone. The 
ionization in the experimental flame was increased by 
potassium seeding which is being carried out by the 
seeding apparatus: shown in @fig. (/). A solutiongor 

KOH was placed in three neck flask. Two electrodes 
carrying a high voltage were placed at the two extreme 
necks - the cathode being placed deep into the solution 
and the anode remaining a few mm above the solution. 
About 15 ky ac was applied at the two electrodes 
through an auto-transformer and a high voltage step-up 
transformer. The electrical spark between the anode 


and KOH solution surface started around 15 kV and once 
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the sparking started the voltage could be reduced to 
about 10 kV just to maintain the Spark, The spark 
atomized the KOH solution and potassium particles were 
carried away by the flowing air to the burner. Due 

to the spark, some water vapour was alsoncarr edapy 
the outgoing air from the spark chamber to a cooling 
System, Water yapour condensed in the cooling chamber 
and dry air with potassium entered the burner. A 
1/16" diameter 2% thoriated tungsten rod was used as 

a cathode and 1/8" diameter 2% thoriated tungsten rod 
as anode. By changing the concentration of solution, 
the level of seeding was altered. To change the con- 
centration, the entire solution was changed and a new 


solution of different concentration poured into the flask, 
3.4 STRUCTURE OF THE LANGMUIR PROBE 


Throughout the experiment, an air cooled brass 
probe with an outer diameter of 1/16" and an inner dia- 
meter of 1/32" is used. High pressure air was passed 
through the probe to keep it cold. The air pressure 


could be regulated from 0 to 15 lb/sq. inch. 


3.5 STRUCTURE OF THE GRID 


The grid structure shown in fig. (8) is made with 
small brass tubes of outer diameter 1/16" and inner dia- 


meter of 1/32" bent to 90° at both ends and inserted 
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into bigger tubes of 1/4" diameter. Soft solder is 
applied at the junction points. One of the sides of 
these bigger tubes is closed and ‘polyflow tubes of 
3/8" diameter is connected to the open ends using a 
Ae D3] 3" junctions. High pressure air of about 


70 lb/sq. inch flows through the grid to keep it cold. 
3.6 ELECTRON DENSITY MEASUREMENT 


Two methods are used to measure the electron den- 
sity of the experimental flame, namely the "saturation 
grid current" and the “negative probe for ion current". 
The grid described in Section 3.5 is placed on the flame 
plasma and is biased negative with respect to the 
burner with the help of a fluke 412-B high voltage power 
supply, the polarity, of which can be changedswithe res 
pects tovground, | The voltagesis varied ands thescurrent 
is measured by a Keithley 153 ammeter until a saturated 
grid current is obtained. At low grid voltage, the 
sheath around the probes remains small and the current 
increases with the increase of voltage. At higher 
yoltage, the sheath formed around the probes of the 
grid expands to such an extent that the sheath overlap 
each other. As the voltage is increased further, the 
sheath does not expand in a horizontal way, but expands 
like a planar sheath in the vertical GLTeCCLOn el lnc 


rate at which the ions are formed is equa. tLoOptnesrace 


“eile acccibe aioe een 
aylawinige” his 1 aink 
ewe'hhs : me Pree 
eT): a pera Pan 
ut} ud BGG ies ageleti ash tans 
ee > Cae tag 
eps at omnia a cant sey 


Seerw;; © ee 


&: quadane ay eo; 2 
ae ee ona ih 4). cerite 
reer ry al oherineh ale. deiv 


49 


at which they are being attracted by the grid. Consequently, 


the grid current saturates with the grid voltage. 


The grid is moved to different parts of the flame 
and the electron density is calculated from the satur- 
ated grid current using the formula I = enveA where 
A= Cross section area of flame, An electron density 


eine is calculated from the saturated grid 


Oren <i) 0 
current which was placed at a height of 10 cm from the 
top of the burner. The current-voltage characteristics 


of the grid from which the electron density is being 


calculated is shown in fig. (9). 


the -eleceron density is also measiredsbyep lacing 
annegativelysbiased probe in the plasma, the electron 
density being calculated from the equation (4) provided 


chigt Raft, is greater than 2. 
3.7 MEASUREMENT OF ION DENSITY IN A PLANAR SHEATH 


Measurements in a thick planar sheath are accomp- 
lished by placing an air cooled grid biased to -400V 
with respect to a burner about 8 cm above the flame. 
The electron density at that height of the flame is 
calculated from the saturation grid current method 
described in Section (3.6) and is found to be v2 x Tdeaetae: 
A sheath of about 2 cm thickness is generated below 
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1 10 100 300 


GRID VOLTAGE (VOLTS) 
The grid current as a function of grid voltage measured at a 


distance of 10 cms from the top of the burner. 
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probe is biased positive with respect to the grid with 
the help of a Fluke power supply and a Keithley 153 
microvolt ammeter is used as a null point detector. 
The probe voltage is then increased until the null point 
is obtained and then the probe is biased -10V with 
respect to this space potential. The current collected 
by the probe at that negative voltage gives the in- 
formation about the ion density at that particular 
point in the sheath. The probe is then moved up and 
down with respect to the Saal cae the same procedure is 
repeated. The ion current is measured inside the plasma 
in the same way. This result is shown in fig. (10). 
iiewms pace, potential “inside. and Oucs ties tie=sheacieds 
measured using the null deflection technique. This 
Homa sorsnownm tn) tg. @CL))s. 
3.8 VELOCITY AND TEMPERATURE MEASUREMENT AT THE TOP OF 
THE CHIMNEY 

The temperature at the top of the chimney where 
the probe was placed is measured a Leeds and Northrup 
Optical pyrometer, no 8621-C, and then compared to 
the thermocouple measurement. A copper-constantine 
thermocouple is used. 

The yelocity V, of the gas is calculated from the 
tracks of illuminated Aluminum oxide which is injected 


at the bottom of the chimney and then photographed with 
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The variation of ton current inside and 
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FIGURE (11) The circuit diagram for the measurement of space 


potential and ionization density in the flame. 
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~ 


a known shutter speed. The shutter speed is measured 
using a light source, photodiode and oscilloscope. 

A mean value of the distances travelled during the 
shutter opening is measured and the VeVOcimey. Ve is 
calculated from the known shutter speed. Thi see Tocrcy. 


is found to be equal to 185 cm/sec. 


aye, 


CHAPTER IV 


' DISCUSSION 


The measurements were carried out in a plasma to 


Reg 
verify that with the appropriate variation of —°- 


the current did indeed change from the Rei cain ao 
current to the field limited current. Measurements 

were also carried out in a thick planar sheath. The 
latter measurements had two purposes, the first to 

check the feasibility of convection measurements in 

the sheath, the second to utilize a useful property 

of the sheath - the ability to produce a large change 
DieLonmvelocity: (i.e. in UE ) while at the same time 
keeping the product (ion density x ion velocity) constant. 


The sheath cofvection current I = F(n.e v.) under these 


a8 
conditions would remain constant while the field limited 


current would vary with the ion density. 


The ion density at different portions of the flame 


was determined from the relationship n_ = = where A = 
fe) ev -A 
cross-sectional area of the flame, i = the total satur- 
ated grid current. The ionization density was observed 
to remain fairly constant throughout the flame except 
in the combustion region of the flame. By working at 
different distances from the burner top, it was possible 
to sample a range of ionization density of from 


oh os 1044 sm? (at ~10 cm from the top of the burner) to 


v2 x 1022 /m> (at ¥30 cm from the burner at the top 
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of the glass chimney). 


Current-voltage characteristics shown in fig. 
(12 a,b) were measured at lower densities with the probe 
Situated above a glass 'chimney!' which confined the hot 
gases. The characteristics were calculated from equation 
(Za)ebyeputting the value of the plasma density Mos cal- 
culated from the saturation grid current when the 
density is 2 x Wire ane and from the ion probe current 
using equation (5) when the density is 2 x Lopes 
The theoretical and experimental current showed fairly 
good agreement at the two densities 2 x WOre ne and 
2aex Ome ne whereas the thick sheath theory of Clements 


(20) 


and Smy and the well known static continuum theory 


(5) 


of Su and Lam show very large errors. 


In the case of probe potentials ranging from 1 to 
10V the experimental points did not follow the theoretical 
curve because in this low probe potential region the 


value of log (r6/tp! remained less than l. 


To yerify the theory of probe measurement in a 
thick sheath condition, a planar sheath of about 2 cm 
was generated below the grid at a density of 2 x 1077 /m?>, 
By considering the flux of ions up the flame to be 
constant n, G1, Et Ve) i, it. is possible to vcalculate 
the electron density ny at any point in the sheath from 


the measured E-field at that point. The E-field was 
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FIGURE (12-a). The, current: yoltage,characteristics of aycy- 
lindrical probe. The solid line shows the theor- 
etical probe.current according eto, equation, (24) 


and the circles are the experimental probe current 
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rial l at iene. cond Ieee tae feud 
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The current voltage characteristics of a cylin- 


drical probe. The solid line shows the theoreti- 


cal probe current according to equation (24) and 
the circles are the experimental probe Curreng.av 
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calculated by determining the Space potential of the 
plasma, at every point inside and outside the Sheath. 
The spacial determination of the space potential at the 
Sheath edge was very difficult, because the sheath 

edge vibrated slowly. Consequently, at every instant 
the potential changed. Again, the determined space 
potential was not very accurate because the + 5n.A. 
scale was used for the null current technique. Le 

was found that the E-field follows the well known space 
charge equation which decreases away from the grid. 
From figs. (10)and(13 )the measured ion current inside 
the sheath was observed to be inversely proportional 

to the electric field. The ion current also increases 
at a constant rate through the sheath but’ saturates 
close to the grid and sheath edge. The reasons for 
thismsaturation 15 not clearly understands 


The measurements of fig. (14) show some of the 


(ae and together 


earlier measurements of Clements and Smy 
WeEenecomeroL chose. of fig. Gl2"a,b) 20 wall cureeneseare 
normalized with respect to the calculated sheath/ 
convection current. These normalized currents are shown 
as a function of the parameter Oe . The theoret- 
peal lyscaleulLated EEN ene ce currents are shown 


as a horizontal line and the field limited current by 
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0.01 


BEGiRE (G4 )aethe ceransrtion sfrom sheath/convection to field limited 


current of a cylindrical probe shown with ther VarTiatroneoL oe 
vx 


@ Clements and Smy from ref-21 
In plasma at the top of the chimney 


& Inside the sheath 


The horizontal solid line is calculated from the sheath/conyection 


theory and the solid sloping line is from field limited current 


theory. 
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the sloping line. At higher densities, 4 x Omen to 
Bh BX 10>" /m-, Clements and Smy observed that the probe 
follows the sheath/convection current theory and their 


experimental points are shown by Ser Cleese et mela) 


; Oe} 
which shows good agreement in the range of S- <1, 
O, ¥X 
When === exceeds 1, the probe current does not 
VX 

follow the sheath/convection theory. At low ionization, 
romecy 
= becomes much greater than unity and the 
7x 


ions entering into the sheath do not reach the probe 
because of the weak E-field present in the sheath, 

The field limited theory shows good agreement with the 
experiments; those performed at low densities are 
indicated by the squares, and those measured inside the 


sheath are indicated by the triangles. 


The ion mobility used in these calculations is 
taken from the most recent value for chemi-ionization 
given by Bradley and Ibharim ce Ory Sex on m*/V-sec 
at atmospheric pressure and 2000° K. The mobility 
is scaled with temperature from the measurements per- 
formed at low temperature. The average flow velocity 
of the plasma is calculated from the statistical mean 
of the particle trace length on the photographic plate 


exposed at a known shutter speed. 


The technique of measuring the variation of ion 


current with probe bias to determine the ion density 
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would seem to be a useful one. 


From the point of view of ionization measurements 
the experimentalist has enough variables at hand to 
avoid the region log Mey a < 1 (namely bias voltage, 
probe. radius, probe/plasma velocity tore al louie nigh 


ion densities or large diameter water cooled electrodes. 
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CHAPTER V 


Conclusions 


From the results shown in TUGS.0, bo) oltre sec loanr 

that for Rea > ih , the sheath/convection formula 
is See ee whereas the field limited current agrees 
reasonably well with the experiment. The sheath measure- 
ments not only show agreement with theory but also 
provide convincing proof that the space charge limited 
current model I = 1(noV¢) is not applicable for a 

c | SX 
Since the slope of the plots would be zero in this case. 
Together with the earlier results of Clements and 


any 


the measurements show that it is always the 
lesser of the sheath/convection and field limited 


currents which is appropriate. 


These measurements, of course, do not constitute 
a comprehensive test of the field limited current 
model. In particular, no measurements have been carried 
out with large probe or spherical probe geometries. 
However, the experimental results given here do provide 


strong evidence for the basic principles of the theory, 


It is interesting to consider some implications of 
these theoretical and experimental results. First, 
it is evident that at low plasma velocities, the sheath 
convection model is adequate down to densities of 
ot /m>. Since there exists substantial theoretical 


and experimental evidence to show that the model can 
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work at densities as high as » Wo a (Clements, Kerr, 


Offenberger and smy) {3) we see that the sheath/con- 


vection model in itself is in fact of broad application. 
At higher velocities, the Roo AGS ine catron 
yx 
becomes more dominant. For instance, in shock tube 
or plasma jet studies near atmospheric pressure, velocities 


) 
~ lO°m/sec or more are common. With flame velocities 


~ 1 m/sec we see that the increase in Ry will be a 


2) 
re 


£actor OF a167 necessitating an increase in no O 
of ioe (tee ves) lige Yee Ne Thus it can be expected 
that the field limited current model be of considerable 
application vat high) velocities. “lias important «to 
note here that the parameter Ryo is not affected by 
the probe dimensions and thus we cannot move into a 
preferred regime of operation by merely changing the 


probe dimensions. 


The equations for cylindrical sand soheried! 
probes in a moving plasma for the situation where the 
probe current is field limited, are deduced. No comp- 
rehensive experimental check of these equations has 
been made. However, it has been established that with 
increasing coils , the transition tor fielaglamired 

x 


current does occur at the Fight, poine eandsecnesCcurrent 


is of the expected magnitude. 
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SUGGESTIONS FOR FUTURE WORK 


The recombination co-efficient of the plasma 
seeded with any other alkali material can be calculated 
using these probe theories to measure the electron or 
ion densities with respect to time after collapsing a 
planar sheath produced in the described way. Again 
the transition from the sheath to the plasma can be 
calculated using the same probe theory at least in the 
case of weakly ionized plasma... The transit time of 
ions can be determined by vibrating this type of 
planar sheath, superimposing an ac voltage ona de 


voltage. 
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